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' We calculate prompt photon production in high-energy nuclear collisions. 

■ We focus on the broadening of the intrinsic transverse momenta of the partons 
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. in the initial state from nuclear effects, and their influence on the prompt 
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photon pt distribution. Comparing to WA98 data from Pb+Pb collisions at 
y/s = IIAA GeV we find evidence for the presence of nuclear broadening at 



r"| ' high Pt in this hard process. Below pt ^ 2.7 GeV the photon distribution is due 

to small momentum transfer processes. At RHIC energy, y/s = 200A GeV, the 
. effect of intrinsic transverse momentum on the spectrum of prompt photons 



is less prominent. The region pt = 3 — 4 GeV would be the most promising for 



. studying the nuclear broadening effects at that energy. Below pt = 2 — 3 GeV 

the contribution from large momentum transfers flattens out, and we expect 
that region to be dominated by soft contributions. 



I. INTRODUCTION 

The WA98 collaboration recently reported data on direct photon production in lead on 
lead collisions at CERN SPS energies (^s = 17 AA GeV) at ~ 1.5 - 4 GeV g. They 
made the interesting observation that in central Pb+Pb collisions the multiplicity of direct 
photons per nudeon-nucleon collision in that range of pt is enhanced relative to proton- 
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proton collisions at similar energy. (For a review of data and calculations of single photon 
production in proton and pion induced scattering off various targets see [0].) It is the 
purpose of this paper to analyze to what extent, and in which range of pt, that data can 
be understood within perturbative QCD (pQCD), after introducing an intrinsic transverse 
momentum of the partons. That transverse momentum is due to the partons being confined 
in the initial-state nucleons, gluon bremsstrahlung, as well as from multiple soft scattering 
of the nucleons prior to the hard scattering. It can be identified in the nuclear dependence 
of the intrinsic transverse momentum of the partons predicted in and observed recently 
for dimuon production in p + A collisions |Q. 

Thus, the purpose of this paper is to analyse the direct photon data of ref. [|l| to probe 
nuclear broadening effects on the transverse momentum distribution of direct photons. Those 
effects are expected to increase the prompt photon cross section strongly, because a part of 
the Pt of the photon can be "supplied" by the incoming partons rather than in the elementary 
semi-hard scattering itself 0. 

Even so, we explicitly restrict the computation to include only the contribution from 
parton-parton scattering amplitudes at a minimum momentum transfer of > 1 GeV. Thus, 
comparing the calculated pt distribution of photons to the data, one can determine the pt 
scale below which production of single photons is dominated by reactions at small momentum 
transfer, < 1 GeV. That is probably the region where a pQCD description looks suspicious. 
Note that without imposing a cut on the minimum momentum transfer in elementary parton- 
parton scatterings the perturbative contribution becomes very large at low pt, precisely 
because of the fact that pt ^ 1 — 2 GeV photons can be produced with almost no momentum 
transfer, the pt being supplied by the intrinsic kt of the partons. Therefore, care must 
be taken when comparing perturbative computations accounting for the violation of the 
DGLAP approximation due to parton intrinsic transverse momentum to data. This primary 
transverse momentum is to large extent due to the gluon bremsstrahlung. With increasing 
energy, primary transverse momentum plays a less important role because a lot of gluons are 
radiated "long before" the hard collision [^]. Thus the distribution of partons near the point 
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of the hard colhsion will be determined mostly by gluon bremstrahlung, which is universal 
- independent of the target. 

We shall also make predictions for Au+Au collisions at nominal BNL-RHIC energy, ^/s = 
200A GeV. Prompt photon spectra at RHIC will be measured in the future by the PHENIX 
collaboration. We shall discuss the importance of intrinsic transverse momentum at that 
energy as well, and identify the pt where the photon distribution from pQCD interactions 
with momentum transfer > 1 GeV flattens out. 

The significant effects from nuclear broadening of the transverse parton distribution on 
photon production in A + A collisions are in general agreement with the description of the 
yield of electrically neutral pions in 0. It proved essential for the understanding of 7r° 
production in Pb+Pb at SPS energy g] {^/s = UAA GeV; pt = 1 - A GeV) to account 
for the pf-broadening. Also, measurements of both the prompt photon and pion production 
cross section in future p + A experiments at BNL-RHIC can determine the scaling of the 
nuclear broadening with A 0,^]; that could provide additional information to that obtained 
from the Drell-Yan process and J/\E', T production. 

II. PROMPT PHOTON Pt DISTRIBUTION IN P + P 

We consider the contributions from Compton-like scattering {g+ V— ^ 7+ V), annihila- 
tion {q+q —>■ 7+5'), plus collinear bremsstrahlung off a (anti-)quark produced at midrapidity. 
Assuming the applicability of the QCD factorization theorem, the corresponding expressions 
for the Compton and annihilation subprocesses inp + p— >7 + X are 

= J dXadx6d^A;tad^/ct6/(fcta)/(fctt)Ga/A(a;a,Q^)Gfc/B(a;6,Q^)^^5(s + t + u)eo • 

(1) 

s, i, u denote the Mandelstam variables for the a + 6 — >■ 7 -(- c elementary process. The sum 
extends over all possible partons in the initial state, i.e. gluons and m, d, s ( ant i-) quarks, 
and over any final state interactions, da/dt denotes the elementary hard-scattering cross 
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section for the corresponding process, averaged over all possible spin and color orientations 
in the initial state, and summed over those in the final state. Explicit expressions can be 
found e.g. in 0. 

The Ga/A{xai Q"^) are the twist-2 infrared dominated matrix elements, for which we em- 
ploy the CTEQ4L parametrization |T^. (We consider only the processes with unpolarized 



particles, so twist-4 would be the next to contribute.) At rather small-x, nuclear shadowing 
of the parton distribution functions may result in considerable suppression of the gluon and 
sea-quark distributions at moderate values of |[Tl|. This may then in turn affect the 
Pt distribution of prompt photons at collider energies (RHIC, LHC) in the pt-range of a 



few GeV 0,|T1. On the other hand, at CERN-SPS energy, ~ 17 - 20A GeV, and for 



~ 1 GeV, we rather deal with the parton distributions at x ~ 0.1, which is in the domain 



of antishadowing effects. The antishadowing effect is rather small, however: EKS [jTj] and 



FGS [jT5[ estimates predict about 10% — 15% enhancement. As we shall see below, this ap- 



pears to be well within the uncertainties of the value of the intrinsic transverse momentum, 
and thus will not be considered here in more detail. 

The functions f{kt) parametrize the transverse momentum of the partons in the initial 
state. For simplicity, we assume a Gaussian distribution. 

The distribution contains one tunable parameter, namely the average intrinsic transverse 
momentum of the partons in the initial state, (k'f). We shall present results for various values 
of (k'f) below. The limit (k^) — * recovers the usual coUinear factorization. 

The intrinsic kt for plane waves bound in a nucleon can be estimated by applying the 
uncertainty principle: 

~ ^ ~ 0.37 GeV . (3) 

Here, r^r ^ 0.85 fm is the radius of a nucleon. Initial state gluon bremsstrahlung can give 
rise to larger intrinsic kt on the order of 1 GeV. 



For the factorization scale we assume = {2pt)'^, and as we restrict ourselves to the 
leading logarithmic approximation, we assume that the renormalization scale is the same. 
We employ the one-loop expression for the running coupling constant as(Q^) for Nf = 4, and 
Aqcd = 236 MeV, as appropriate for CTEQ4L. Collinear divergencies arising from partonic 
processes with large intrinsic transverse momentum are cut off by 

00 = 6162 e _ iGeV^) , (4) 

where we define 



" p +T+ • 



In this way, only subprocesses with momentum transfer > 1 GeV are taken into account, 
for which one can hope that factorization is a reasonable assumption. The functions Bi and 
G2 ensure that parton a is moving to the right, while parton h is moving to the left, 

81 = Q{Ka) , 62 = . (6) 

Using on-shell kinematics, the four-momenta of the incoming partons are given by 

/l P - 1 k"^ \ 

(i^.4,M= (^x.v/^+^,4,-i..v/^+^) . (8) 

After using the (5-function in eq. (|I|) to solve for x^, the remaining integrals can be performed 
numerically via Monte-Carlo techniques. Numerical convergence can be checked by compar- 
ing spectra obtained with varying sequences of random numbers and increasing the number 
of Monte-Carlo points. 

Collinear bremsstrahlung is computed as follows. First, the cross-section for (anti-)quark 
production at Ucm = with transverse momentum ptc is computed in the standard fashion, 
as described above. This is then convoluted with the QED collinear fragmentation function 
ZcD.ilq{zc) = a[l + (1 - ZcY] log(p^/AQCD)/2vr into 7 + X, where Zc = Pt/ptc, 
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-E"^^ = J dxadxbdzcd^ktad'^ktb 

s dcr 



f{kta)f{ktl>)Ga/A{Xa, Q')G'vb(x6, Q')D^/,{z,) — -j6is + t + u)Qo . (9) 

TCZ^ dt 



Here, we employ = {2ptc)'^ for the factorization and renormalization scale of the hard 
process. For the bremsstrahlung contribution we do not sum explicitly over all hard sub- 
processes leading to the production of a (anti-)quark at midrapidity but rather employ the 
method of effective structure function, 

where ej is the fractional electric charge of the splitting (anti-)quark. Further, da/dt = 
TT{as/sy{s'^ + u^)/i^. The 0-function is defined as in (§J^). Note that in these expressions 
s, t, u refer to the partonic subprocess, "before" splitting of the final-state (anti-)quark into 

It is common practice to multiply the r.h.s. of ([I|) and by a i^-factor in order to 
account for NLO contributions. Full NLO computations have been performed []TB[ for the 
case without intrinsic kt of the partons and indicate ~ 1.8 — 2.0 for p + p collisions at 
energies around ^/s = 20 GeV, and for photon transverse momenta in the range 2 — 4 GeV. 
Our major interest is in the relative yields of prompt photons in p+p and in A + A collisions, 
and in that sense the i^-f actor plays no essential role. In any case, we shall fix the value of 
K from p + p and then keep its value for A + A. 
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FIG. 1. The differential direct photon cross section vs. pt (at center-of-mass rapidity) for pro- 
ton-proton collsions at ^/s = 19.4 GeV. 



In Fig. ^ we compare our results for p + p ai ^/s = 19.4 GeV, for various intrinsic (k^) 
and = 2, to E704 data fl^. Clearly, the calculation assuming pure DGLAP evolution, 
i.e. (k^) = 0, falls short by almost an order of magnitude. The best fit is obtained using 
(kf) = 1.3 — 1.5 GeV^. At first sight this appears somewhat larger than the "standard" 
(kf) = 0.9 GeV^ from |]5|,p!8||. However, recall that we employ a larger = 4pf scale 
than P,[18|, who used = p1/2. In accord with the picture that the intrinsic transverse 



momentum is due mainly to gluon radiation in the initial state, a larger scale therefore 
has to correspond to larger {k"^). This general expectation is confirmed by our results. (The 
value for (A;^) therefore only acquires physical meaning together with a definition of Q"^.) 

The above calculation serves to fix the i^-factor and the strength of the DGLAP violation 
on the p + p level. From there we conclude that with K = 2, = Ap^, a value of 
{k^) ^ 1.3 — 1.5 GeV^ is required to fit the data at that energy. 
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III. HARD PHOTON PRODUCTION IN NUCLEUS-NUCLEUS COLLISIONS 



In the previous section we discussed hard photon production m p + p or N + N colhsions 
(we neglect the differences between protons and nucleons in our calculations). To extend 
our calculations to nucleus-nucleus collisions we have to substitute the parton distribution 
functions Gn{x, Q) in eq. (|l]) by Ga{x, Q, b') for the projectile and Gb{x, Q,b — b') for the 
target nucleus. Here, Ga{x, Q, b') = A ■ TA{b')GN{x, Q) and Gb{x, Q,b-b') = B - Tsib- 
b')Gi\f{x, Q) [|l^], with b the nucleus-nucleus impact parameter, and b' the distance from the 
center of the projectile. Note that b and b' are two-dimensional vectors. 

The nuclear thickness functions TA^sip') are defined as 

TA,B{b') = J dzpAAb',^) with J d'b'dzpib',z) = 1 . (11) 

A (B) are the number of nucleons in the projectile (target). We account for the dependence 
on the nucleus nucleus impact parameter b as due to the geometry of the colhsion only. If 
one integrates over all impact parameters to compare with minimum bias data, this yields 

<^ = AB ■ <^ . (12) 

To compare to the experimental data by the WA98 collaboration we need to calculate the 
direct photon cross section for their centrality trigger. To take into account the dependence 
of the nuclear parton distribution functions on the geometry we employ the semiclassical 
approximation of the Glauber-model that was developed for hadron nucleus collisions in 



ref. [BOl and extented in ref. ETI to nucleus nucleus collisions. 



The inclusive spectrum of hard photons in A + B collisions at the impact parameter b is 

Ey and are the energy and the momentum of the emitted photon, and is the inclusive 
cross section for direct photon production. 

Tab is the nucleus-nucleus thickness function. For large nuclei with a slowly varying 
density distribution Tab is given by 
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TAB{b) = J d'b'TA{b')TB{b - b') . (14) 

In a. A + B collision at impact parameter b, the nucleus-nucleus thickness function allows 
only collisions of partons in overlapping slabs, i.e. of partons with the same coordinates in 
the plane transverse to the beam direction. 

Finally, the differential number of direct photons in nucleus-nucleus collisions integrated 
over impact parameters < 6 < 6^, is given by the ratio of the differential cross section to 
the inelastic nucleus- nucleus cross section, both being integrated over the region [0,6™.]: 

diV^ _ E, J^- d'b'^jb) _ E.^J^-d^bAB-TMb) 



From this equation, we can also obtain the prompt photon multiplicity for p + A scattering. 
To that end, the nucleus-nucleus thickness function T^Bib) has to be replaced by the nucleus 
thickness function TA{b), and B = 1 throughout the equation. For p + A we have bm = oo, 
corresponding to minimum bias scattering. 

In our calculations we used the Woods-Saxon density distribution of ref. and assumed 
= 30 mb at y/s ~ 20A GeV. That yields a total inelastic cross section for Pb+Pb of 
(jfJl^^ = 7100 mb. crfJt^^ is given by the denominator of eq. ([T5|) with km = oo. For large 
values of AB one can approximate the denominator by 

/ d'b (1 - (1 - T^B(6)aD^'') ^ / d'b (1 - exp(-AET^B(6)ar )) • (16) 



In ref. |21[] this is called optical limit for the inelastic cross section. 

The geometrical factor in eq. (|1^) which multiplies the p + p differential cross section to 
yield the photon multiplicity equals 24.3/mb for the 10% most central Pb+Pb collisions, as 
analysed by WA98 (corresponding to bm ~ 4.8 fm in our calculation). For minimum bias 
p+C it is 0.056/mb, and 0.127/mb for p+Ph. These values allow to calculate the yield of 
prompt photons in p + A collisions at SPS energy from the cross section in p + p as shown 
in Fig. 0. 

For nominal BNL-RHIC energy, ^/s = 200A GeV, using a^^'^ = 40 mb, we obtain for 
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that factor 0.119/mb for p+Au, 0.048/mb for p+C, 0.062/mb for p+S, and finally 23.1/mb 
for Au+Au (with 10% centrahty cut; af^^^"" = 7024 mb; ^ 4.7 fm). 

IV. NUCLEAR BROADENING OF THE INTRINSIC Kt DISTRIBUTION 

The nuclear broadening of the intrinsic transverse momentum distributions of the partons 
was observed experimentally by the broadening of the transverse momentum distribution of 
final state particles. For example, ref. ||2^ analyzed pt distributions of mesons made of light 
quarks (vr, K, . . .); ref. analyzed Drell-Yan pairs and quarkonium states {J/ip and T). 

The aim of this paper is to discuss whether that broadening can explain the discrepancy 
of the WA98 data for Pb+Pb collisions to simple extrapolations of proton-proton data and 
calculations. We define the broadening as 

Ak^ = {k^)AB - {k^)pp , (17) 

where {kf)AB(pp) is the average transverse momentum in A + B {p + p) collisions. For 
Pb+Pb collisions values of Akf = — 1 GeV^ will be used in our calculations below. Larger 
broadening, as suggested for example by ref. |^ for hot QCD matter, seem to be ruled out 
at SPS energies as will be shown in section |V|. The value suggested in for cold (hot) 
QCD matter is 

Akf = 0.2 (3.0) GeV' ■ — ^ (18) 
* ^ ^ lOfm ^ ^ 

L is the length of the QCD medium through which the parton propagates before the hard 
process. The average thickness of a sphere is the volume divided by the area: 

{L),phere = {^71 R^) / (tI R^) = ^Ra,B • (19) 

For proton-nucleus collisions that value has to be divided by two, since we are interested in 
initial state interaction only (the outgoing photon does not interact strongly). For nucleus- 
nucleus collisions another factor two is necessary, since both the intrinsic transverse mo- 
mentum of the projectile and of the target can increase due to soft inital-state interac- 
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tions. Ra,b denotes the radius of the projectile (target) nucleus. For Pb+Pb this means 
L ^ 2 ■ |i? = I ■ 6.6 fm = 8.8 fm. From eq. ([T8|) we thus obtain for cold (hot) QCD matter 

Ak'^ = 0.2 (3.0) GeV^ ■ 0.88 = 0.176 (2.64) GeVl (20) 

For gluons Ak^ should be larger by the ratio of the Casimir operators of the octet and triplet 
representations of color-SU(3), = 9/4. While these number are only rough estimates of the 
order of magnitude, we shall see below that Akf ~ 0.5 — 1 GeV^ is in the range needed to 
explain the WA98 prompt photon excess for ~ 2.5 GeV. 

Such small broadening is also compatible with the fits to p + A data at Fermilab ener- 
gies 0]. Those fits yield 

Ak^ = 0.133 (0.027) GeV^ ^(^^) ' - 1 j (21) 

for J/ip (Drell-Yan pair) production. For p+Ph this is Ak^ = 0.49 (0.10) GeV^. Assuming 
linear scaling |2^ of Ak^ with the length traversed through the QCD matter, L, the above 



value has to be multiplied by two, for minimum bias lead on lead collisions: 

AA;i2(Pb + Pb) = 0.98 (0.2) GeVl (22) 

Regarding the Drell-Yan process, the Akf broadening of the quarks in the initial state equals 
that of the muon pair in the final state. In our case, however, the dominant processes involve 
a gluon in the initial state, and so the broadening in the initial state is expected to be larger 
by a factor (1 + 9/4)/(l + 1), i.e. 0.2 GeV^ 0.325 GeV^. For a central trigger, L and thus 
Akf could even be somewhat larger. 

To summarize, values for Akf up to 1 GeV^ appear to be within the range given by 
various fits to data, and shall be employed below for the computation of prompt photon 
production in nuclear collisions. 
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V. RESULTS FOR CERN-SPS ENERGY 




Pt [GeV] 

FIG. 2. The transverse momentum distribution of prompt photons at center-of-mass rapidity 
in the 10% most central Pb+Pb reactions at ^/s = 17 AA GeV. The data is compared to 
calculations for various intrinsic (/cf). 

Fig. 1^ depicts the differential multiplicity of direct photons vs. transverse momentum in 
central Pb+Pb reactions at ^/s = 17 AA GeV. We compare calculations with various values 
for the intrinsic transverse momentum, (kf), all corresponding to = 2. The data are from 
ref. Again, as was the case f oi p + p collisions at similar energy, the calculation assuming 
pure DGLAP evolution underestimates the data strongly. However, the "standard" intrinsic 
transverse parton momentum of (k^) ^ 1.3 GeV^ extracted already on the p + p level (see 
section |IlD improves the agreement a lot. Within the experimental error bars, the WA98 data 
leaves room for up to Akf ~ 1 GeV^ additional broadening from nuclear effects. However, 
we can also conclude that transverse momentum broadening much beyond Akf ~ 1 GeV^ 
is not seen in the WA98 data. 

Another interesting result is that the calculated photon spectrum below pt < 2.5 GeV 
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depends only weakly on the amount of intrinsic transverse momentum. That is because 
of the cut-off on the momentum transfer, eq. (^). Therefore, the excess photons seen in 
the data, as compared to our calculation, can be attributed to small-momentum transfer 
scatterings (final-state interactions). The computation of that contribution to the photon 
yield is, however, out of the scope of the present manuscript; see e.g. for recent work, 
and references therein. 

VI. RESULTS FOR BNL-RHIC ENERGY 

We now turn to nuclear collisions at collider energy. Here, we focus on Au+Au collisions 
at ^/s = 200 A GeV and central rapidity, y ^ 0. (Estimates for p + A collisions can be 
obtained from our Au+Au results by an appropriate scaling as explained in section |T|.) At 
relatively small transverse momenta nuclear shadowing effects may reduce the prompt 
photon cross section somewhat ||T2| , p!3| . However, for ~ 2 GeV the effect is not very 
large, < 20%, and we shall not investigate it in more detail here. On the other hand, it has 
been suggested that at high energies the gluon distribution in nuclei may saturate at small 
fractional light-cone momentum x (but much larger x than for protons). This must 
happen at least when the unitarity bound for the inelastic scattering of a QCD dipole off a 
nucleus is reached in the x — plane perhaps earlier. Simply speaking, the nucleus 



becomes a "black disc" in the infinite momentum frame. If heavy nuclei but not protons 
become black, or almost black, jets with pt < Quack will be absorbed in the final state in 
y4 + y4 scattering^, but not in p + p scattering. So high pt jets will be expected for pt > Quack 
only. Thus, qualitatively, the pt distribution should be rather different from p + p at not too 
large pt- The yield of jets will be suppressed but the transverse momentum distribution will 
be wider. 



^They might produce a "saturated plasma" |27| at central rapidity which would then also con- 



tribute to photon production |25 
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For A = 200 nuclei, the unitarity limit for the inelastic gg-nucleus cross section is reached 
at a; ~ 10^^ for Q ~ 1.5 GeV, and at x ~ 10^'^ for Q ~ 2.5 GeV (if nuclear shadowing is 
weak) [O. If shadowing is also taken into account, the a;-values where the unitarity limit 



is approached become even smaller |]I5[. Therefore, at rapidity y ~ and at RHIC energy, 
^/s = 200A GeV, we expect only small effects on prompt photon production. 

At high energy, the bremsstrahlung contribution to the prompt photon yield is sub- 
stantial ||12| ,[T3| , |T6| . In this process, a (anti-)quark is produced first in the central rapidity 
region with transverse momentum ptc (see section |T|). That quark then radiates an almost 
collinear photon with transverse momentum pt = zptc- In principle, the quark could in fact 
suffer further collisions which in turn may suppress the contribution from soft (i.e., small z) 
bremsstrahlung to prompt photon production |[1^J2^ . This is due to the well known Landau- 
Pomeranchuck-Migdal effect. The quark could also suffer some energy loss prior to emitting 
the photon. The radiative energy loss of the quark is mainly due to radiation of soft gluons. 
Therefore, the formation time (or the formation length) of the radiated gluons cannot be 
much smaller than that of the photons. Thus, quark energy loss might not interfere with 
the radition of the photons in the final state. In any case, the complete extinction of all 
bremsstrahlung would represent an upper limit for the sum of both effects . This would 
change our result by less than a factor of two. Therefore, we do not discuss these effects in 
more detail here (see ll^), since their largest possible impact on the photon production cross 
section is smaller than that from nuclear broadening of the intrinsic transverse momentum. 
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<k, >=0. GeV" 
<kt^>=0.9 GeV^ 
<k,^>=1.5 GeV^ 
<k,^>=2.4 GeV^ 




1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 
Pt 



FIG. 3. The transverse momentum distribution of prompt photons at center-of-mass rapidity in 
the 10% most central Au+Au reactions at ^/s = 200A GeV for various intrinsic (/cf). 



In Fig. 1^ we show prompt photon spectra at RHIC energy, for the 10% most central 
Au+Au events at ^/s = 200A GeV. As before, we assumed K = 2. One observes that the 
effect of intrinsic transverse momentum in the initial state is much less prominent than at 
lower energy, in agreement with the results of for pion production. The reason for this 
behavior is that the photon spectrum becomes "harder" at higher energy, i.e. it decreases less 
steeply with pt. Nevertheless, at ~ 3 — 4 GeV the intrinsic parton transverse momentum 
can increase the prompt photon multiplicity by up to a factor of 3. That would be the most 
promising kinematical domain for an experimental study of the nuclear broadening of the 
intrinsic transverse momentum distributions. The region of smaller pt might be dominated 
by soft contributions, as indicated by the flattening of the spectra in Fig. ^. On the other 
hand, at larger pt intrinsic transverse momentum becomes less effective. 
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VII. SUMMARY AND CONCLUSION 



Prompt photon data from p + p scattering at y/s = 19.4 GeV obtained by the E704 
collaboration [1^| reveals that standard DGLAP evolution without intrinsic transverse mo- 
mentum, and for our choice for the factorization scale Q = 2pt, underestimates that data. 
Including an intrinsic transverse momentum of the initial-state partons of ~ 1.3 GeV^ allows 
for a fair description of the E704 data. 

Further, we presented results for prompt photon pt distributions in nuclear collisions, 
assuming A-scaling as suggested by the geometry of the Glauber incoherent approach. At 
SPS energy, a good description of the data for pt ~ 2.7 GeV requires additional intrinsic kt of 
the partons in the initial state, i.e. nuclear broadening of the transverse parton distributions 
by Afcj ~ 0.5 — 1 GeV^ (for the same factorization scale as above). Thus, the WA98 data 
provide evidence for the nuclear broadening effect in the pt distribution of prompt photons. 
Photons appear more appropriate than hadrons to test nuclear broadening of the intrinsic 
transverse momentum because they do not suffer from rescattering in the final state. 

Below Pt ~ 2.7 GeV, both the yield and the slope of the WA98 photon data is not 
described as being due to hard scatterings with momentum transfer > 1 GeV. At ~ 2 GeV, 
prompt photons from hard collisions underestimate the WA98 data by almost an order of 
magnitude, for any amount of intrinsic {k^). Thus, that kinematic domain is dominated by 
small momentum transfer processes. 

At RHIC energy, the effect of the intrinsic transverse momentum is much less prominent. 
Only at ~ 3 — 4 GeV can intrinsic parton transverse momentum increase the prompt pho- 
ton multiplicity by up to a factor of 3. That would thus be the most promising kinematical 
domain for an experimental study of nuclear broadening of the intrinsic transverse momen- 
tum at RHIC. Again, at < 2 — 3 GeV the realistic prompt photon distribution (inch 
intrinsic kt) from large momentum transfer scattering flattens out, that region being most 
likely dominated by soft contributions. 
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